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Abstract 

We review properties of cosmological theories for the variation of the 
fine structure 'constant'. We highhght some general features of the cosmo- 
logical models that exist in these theories with reference to recent quasar 
data that are consistent with time- variation in the fine structure 'constant' 
since a redshift of 3.5. 



1 Introduction 

There are several reasons why the possibility of varying constants should be 
taken seriously [1]. First, we know that the best candidates for unification of 
the forces of nature in a quantum gravitational environment only seem to exist 
in finite form if there are many more dimensions of space than the three that we 
are familiar with. This means that the true constants of nature are defined in 
higher dimensions and the three-dimensional shadows we observe are no longer 
fundamental and need not be constant. Any slow change in the scale of the extra 
dimensions would be revealed by measurable changes in our three-dimensional 
'constants'. Second, we appreciate that some apparent constant might be de- 
termined partially or completely by spontaneous symmetry-breaking processes 
in the very early universe. This introduces an irreducibly random element into 
the values of those constants. They may be different in different parts of the 
universe. The most dramatic manifestation of this process is provided by the 
chaotic and eternal inflationary universe scenarios where symmetries determin- 
ing both the number and the strength of forces in the universe at low energy 
can break differently in different regions. Third, any outcome of a theory of 
quantum gravity will be intrinsically probabilistic. It is often imagined that 
the probability distributions for observables will be very sharply peaked but 
this may not be the case for all possibilities. Fourth, a non-uniqueness of the 
vacuum state for the universe would allow other numerical combinations of the 
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constants to have occurred in different places. String theory indicates that there 
is a huge 'landscape' (> 10^°") of possible vacuum states that the universe can 
find itself residing in as it expand and cools. Each will have different constants 
and associated forces and symmetries. It is sobering to remember that at present 
we have no idea why any of the constants of Nature take the numerical values 
they do and we have never successfully predicted the value of any dimensionless 
constant in advance of its measurement. However, the last reason to consider 
varying constants is currently the most compelling. For the first time there is 
a body of detailed astronomical evidence for the time variation of a traditional 
constant. The observational programme of Webb et al [U [3] has completed 
detailed analyses of three separate quasar absorption line data sets taken at 
Keck and finds persistent evidence consistent with the fine structure constant, 
a, having been smaller in the past, at z = 1 — 3.5. The shift in the value of a 
for all the data sets is given provisionally by /S.a/a — (— 0.57± 0.10) x 10~^. 
This result is currently the subject of detailed analysis and reanalysis by the 
observers in order to search for possible systematic biases in the astrophysical 
environment or in the laboratory determinations of the spectral lines. So far 
it has not been undermined or confirmed by other observations (for the most 
reason discussion of the status of uncertainties, see [1]). 

The first investigations of time-varying constants were those made by Lord 
Kelvin and others interested in possible time-variation of the speed of light at 
the end of the nineteenth century. In 1935 Milne devised a theory of gravity, 
of a form that we would now term 'bimetric', in which there were two times 
- one {t) for atomic phenomena, one (r) for gravitational phenomena - linked 
by T = log(t/to)- Milne [5] required that the 'mass of the universe' (what we 
would now call the mass inside the particle horizon M « c^G~^t) be constant. 
This required G cc t. Interestingly, in 1937 the biologist J.B.S. Haldane took a 
strong interest in this theory and wrote several papers [6] exploring its conse- 
quences for the evolution of life. The argued that biochemical activation energies 
might appear constant on the t timescale yet increase on the r timescale, giving 
rise to a non- uniformity in the evolutionary process. Also at this time there 
was widespread familiarity with the mysterious 'large numbers' 0(10'*'') and 
O(10®°) through the work of Eddington (although they had first been noticed 
by Weyl [7] - see ref. [8] and [1] for the history). These two ingredients were 
merged by Dirac in 1937 in a famous development (supposedly written on his 
honeymoon) that proposed that these large numbers (10''") were actually equal, 
up to small dimensionless factors. Thus, if we form N ^ c^t/Gnin ^ 10^°, 
the number of nucleons in the visible universe, and equate it to the square of 
Ni ^ e'^/Grn^ ^ 10*", the ratio of the electrostatic and gravitational forces 
between two protons then we are led to conclude that one of the constants, 
e,G,c,h,mn must vary with time. Dirac [3] chose G oc to carry the time 
variation. Unfortunately, this hypothesis did not survive very long. Edward 
Teller [TU] pointed out that such a steep increase in G to the past led to huge 
increases in the Earth's surface temperature in the past. The luminosity of 
the sun varies as i a and the radius of the Earth's orbit as i? cc G^^ so 
the Earth's surface temperature varies as {L/R^Y^'^ oc G^/* a t^^/* and 
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would exceed the boiling point of water in the pre-Cambrian era. Life would 
be eliminated. Gamow subsequently suggested that the time variation needed 
to reconcile the large number coincidences be carried by e rather than G, but 
again this strong variation was soon shown to be in conflict with geophysical and 
radioactive decay data. This chapter was brought to an end by Dicke [TT] who 
pointed out that the N ^ Nf large number coincidence was just the statement 
that t, the present age of the universe when our observations are being made, is 
of order the main-sequence stellar hfetime, tms ~ {Gm'^/hc)~^h/mnC^ ^ 10^° 
yrs, and therefore inevitable for observers made out of chemical elements heav- 
ier than hydrogen and helium. Dirac never accepted this anthropic explanation 
for the large number coincidences (believing that 'observers' would be present 
in the universe long after the stars had died) but curiously can be found mak- 
ing exactly the same type of anthropic argument to defend his own varying G 
theory by highly improbable arguments (that the Sun accretes material peri- 
odically during its orbit of the galaxy and this extra material cancels out the 
effects of overheating in the past) in unpublished correspondence with Gamow 
in 1967 (see [T] for fuller details). Dirac's biographer has revealed that in 1993 
he expressed 'an article of faith... that the human race will continue to live for 
ever and will develop and progress without limit' [H]. This belief motivates his 
comments relating to the anthropic argument. 

Dirac's proposal acted as a stimulus to theorists, like Jordan, Brans and 
Dicke |13| . to develop rigorous theories which included the time variation of 
G self-consistently, by modelling it as arising from the space-time variation of 
some scalar field (/)(x, t) whose motion both conserved energy and momentum 
and created its own gravitational field variations. The possibility that a varies 
in time has led to the first extensive exploration of simple self-consistent theories 
in which a changes occur through the dynamics of some scalar field. 

2 A Simple Varying- Alpha Theory 

We consider some of the cosmological consequences of a simple theory of time 
varying a. Such a theory was first formulated by Bekenstein [14j as a generalisa- 
tion of Maxwell's equations but ignoring the consequences for the gravitational 
field equations. We completed this theory [15] to include the coupling to the 
gravitational sector and analysed its general cosmological consequences and it is 
referred to as the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) theory below. 
Extensions to include the weak interaction via a generalised Weinberg-Salam 
theory have also been explored, [TBI [T7] .Extensions to include the weak interac- 
tion via a generalised Weinberg-Salam theory have also been explored, [l6l flT] . 

The idea that the charge on the electron, or the fine structure constant, 
might vary in cosmological time was proposed in 1948 by Teller, [10], who 
suggested that a cx (Int)"-^ was implied by Dirac's proposal that G oc 
and the numerical coincidence that a^^ ^ ln(/ic/Gmp), where nip is the proton 
mass. Later, in 1967, Gamow [18] suggested a cx t as an alternative to Dirac's 
time- variation of the gravitation constant, G, as a solution of the large numbers 
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coincidences problem; and in 1963 Stanyukovich had also considered varying a 
in this context [19] . However, any such power-law variation in the recent past 
was soon ruled out by other geological evidence [50] ■ 

There are a number of possible theories allowing for the variation of the fine 
structure constant, a. In the simplest cases one takes c and h to be constants 
and attributes variations in a to changes in e or the permittivity of free space 
(see [21] for a discussion of the meaning of this choice). Thus eo — > e = eoe(x'^), 
where e is a dimensionless scalar field and ep is a constant denoting the present 
value of e. This implies that some well established assumptions, like charge con- 
servation, must give way [22j . Nevertheless, local gauge invariance and causality 
are maintained. 

Since e is the electromagnetic coupling, the e field couples to the gauge field 
as eAf^ in the Lagrangian and the gauge transformation which leaves the action 
invariant is eA^ — > eA^ -|- x,fj.j rather than the usual Af^ + x.^- The 

gauge-invariant electromagnetic field tensor is therefore 

F^, = i iieA,)^^ ~ (eA^),,) , (1) 

which reduces to the usual form when e is constant. The electromagnetic part 
of the action is still 

S,„,^- J d^x^F'^-'F^,. (2) 
and the dynamics of the e field are controlled by the kinetic term 

' " ~2P J (3) 

as in dilaton theories. Here, I is the characteristic length scale of the theory, 
introduced for dimensional reasons. This constant length scale gives the scale 
down to which the electric field around a point charge is accurately Coulombic. 
The corresponding energy scale, hc/l, has to lie between a few tens of MeV and 
Planck scale, ~ lO^^GeV to avoid conflict with experiment. 
The field equations are 

G^. = SttG (T--"- + T^t + T^^^e-'^) . (4) 

The stress tensor of the ip field is derived from the lagrangian £^ — —^d^ipd^ip 
and the ip field obeys the equation of motion 

□ ^ = -e-^'^rem (5) 
CO 

where we have defined the coupling constant oj — {c)/P. This constant is of 
order ^ 1 if, as in [T^, the energy scale is similar to Planck scale. It is clear 
that Com vanishes for a sea of pure radiation since then £em — {E"^ — _B^)/2 = 
0. We therefore expect the variation in a to be driven by electrostatic and 
magnetostatic energy-components rather than electromagnetic radiation. 
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In order to make quantitative predictions we need to know how much of the 
non-relativistic matter contributes to the RHS of Eqn. ([5]). This is parametrised 
by C = LemI P, where p is the energy density, and for baryonic matter £em = 
E"^ j1. For protons and neutrons Cp and C^^ can be estimated from the electro- 
magnetic corrections to the nucleon mass, 0.63 MeV and —0.13 MeV, respec- 
tively [231 . This correction contains the E"^ 12 contribution (always positive), 
but also terms of the form j^a^ (where is the quarks' current) and so cannot 
be used directly. Hence we take a guiding value Cp ~ Cn 10"'''. Furthermore 
the cosmological value of C (denoted C^) has to be weighted by the fraction 
of matter that is non-baryonic. Hence, Cm depends strongly on the nature of 
the dark matter and can take both positive and negative values depending on 
which of Coulomb-energy or magnetostatic energy dominates the dark matter 
of the Universe. It could be that Ccdm ~ ^ 1 (as for superconducting cosmic 
strings, with £cin ~ —5^/2), or Ccdm ^ 1 (neutrinos). BBN predicts an ap- 
proximate value for the baryon density of fie ~ 0.03 (where Ob is the density 
of matter in units of the critical density 3i/^/87rG) with a Hubble parameter of 

= 60 Kms"^ Mpc"^, implying r^cDM ~ 0.3. Thus, depending on the nature 
of the dark matter. Cm can be virtually any number between —1 and -1-1. The 
uncertainties in the underlying quark physics and especially the constituents of 
the dark matter make it difhcult to impose more certain bounds on Cm- 

2.1 The cosmological equations 

Assuming a homogeneous and isotropic Friedmann metric with expansion scale 
factor a(t) and curvature parameter k in eqn. ([4]), we obtain the field equations 
(c^l) 

/ 8ttG / , r „ ,iN r „ ,1 •,2\ k A . , 

[a) = — ('°m(l + CmCxp[-2^]) + p,cxp[-2V'] + -^ J -^ + - (6) 

where A is the cosmological constant. The scalar field obeys 

V:> + 3i/7/' = --exp[-2V;]C„,Pm, (7) 
where H = d/a is the Hubble rate. We can rewrite this as 

(V^a^j = iVexp[-2V'] (8) 

where iV is a positive constant defined by = —2(^^p^^a^ /uj.Note that the 
sign of the evolution of tp is dependent on the sign of Cm- Since the obser- 
vational data is consistent with a smaller value of a in the past, we will in 
this paper confine our study to negative values of Cmi li'^s with our recent 
discussion in Refs. [15]. The conservation equations for the non-interacting 
radiation and matter densities give cx and p^. e"^"^ oc a~^, respectively. 
This theory enables the cosmological consequences of varying e, to be analysed 
self-consistently rather than by changing the constant value of e in the standard 
theory to another constant value, as in the original proposals made in response 
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to the large numbers coincidences. We shall consider the form of the solutions 
to these equations when the universe is successively dominated by the kinetic 
energy of the scalar field ip, pressure-free matter, radiation, negative spatial cur- 
vature, and positive cosmological constant. Our analytic expressions are checked 
by numerical solutions of (O and ([7]). There are a number of conclusions that 
can be drawn from the study of the simple BSBM models with < 0. These 
models give a good fit to the varying a imphed by the QSO data of refs. [2 [3]. 
There is just a single parameter to fit and this is given by the choice 

_ Cn ^ (2±1) X 10-* (9) 

The simple solutions predict a slow (logarithmic) time increase during the 
dust era of fc = Friedmann universes. The cosmological constant turns off 
the time- variation of a at the redshift when the universe begins to accelerate 
{z ~ 0.7) and so there is no conflict between the a variation seen in quasars at 
z ^ 1 — 3.5 and the limits on possible variation of a deduced from the operation 
of the Oklo natural reactor [25 (even assuming that the cosmological variation 
applies unchanged to the terrestrial environment). The reactor operated 1.8 
billion years ago at a redshift of only z ~ 0.1 when no significant variations 
were occurring in a. The slow logarithmic increase in a also means that we 
would not expect to have seen any effect yet in the anisotropy of the microwave 
backgrounds [531125]: the value of a at the last scattering redshift, z = 1000, is 
only 0.005% lower than its value today. Similarly, the essentially constant evo- 
lution of a predicted during the radiation era leads us to expect no measurable 
effects on the products of Big Bang nucleosynthesis (BBN) because a was only 
0.007% smaller at BBN than it is today. 

Theories in which a varies will in general lead to violations of the weak equiv- 
alence principle (WEP) . This is because the a variation is carried by a field like 
ijj and this couples differently to different nuclei because they contain different 
numbers of electrically charged particles (protons). The theory discussed here 
has the interesting consequence of leading to a relative acceleration of order 
10"^'^ [57] if the free coupling parameter is fixed to the value given in eq. ^ 
by using a best fit of the theories cosmological model to the QSO observations 
of refs. [21 [3] . Other predictions of such WEP violations have also been made 
in refs. [28l |23l [29l [30] . The observational upper bound on this parameter from 
direct experiment is just an order of magnitude larger, at 10"^^, and limits from 
the motion of the Moon are of similar order, but space-based tests planned for 
the STEP mission [31] are expected to achieve a sensitivity of order 10~^^ and 
will provide a completely independent check on theories of time- varying e and 
a. This is an exciting prospect for the future. 

2.2 The nature of the Friedmann solutions 

Let us present the predicted cosmological evolution of a in the BSBM theory, 
that we summarised above, in a little more detail. During the radiation era 
the expansion scale factor of the universe increases as a{t) ^ <^/^ and a is 
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essentially constant in universes with an entropy per baryon and present value 
of a like our own. It increases in the dust era, where a(t) ~ t^l^ . The increase 
in a however, is very slow and a ~ 2iVlog(t/ii). This slow increase continues 
until the expansion becomes dominated by negative curvature, ait) ~ or 
by a cosmological vacuum energy, a(t) ~ exp[Ai/3]. Thereafter a asymptotes 
rapidly to a constant. If we set the cosmological constant equal to zero and 
fc = then, during the dust era, a would continue to increase indefinitely. The 
effect of the expansion is very significant at all times. Non-zero curvature or a 
cosmological constant stops the increase in the value of a that occurs during the 
dust-dominated era. Hence, if the spatial curvature and A are both too small 
it is possible for the fine structure constant to grow too large for biologically 
important atoms and nuclei to exist in the universe. There will be a time 
in the future when a reaches too large a value for life to emerge or persist. 
The closer a universe is to flatness or the closer A is to zero so the longer the 
monotonic increase in a will continue, and the more likely it becomes that life 
will be extinguished. Conversely, a non-zero positive A or a non-zero negative 
curvature will stop the increase of a earlier and allow life to persist for longer. 
If life can survive into the curvature or A-dominated phases of the universe's 
history then it will not be threatened by the steady cosmological increase in a 
unless the universe collapses back to high density. 

There have been several studies, following Carter, [52] of the need for life- 
supporting universes to expand close to the 'flat' Einstein de Sitter trajectory 
for long periods of time. This ensures that the universe cannot collapse back 
to high density before galaxies, stars, and biochemical elements can form by 
gravitational instability, or expand too fast for stars and galaxies to form by 
gravitational instability |33l f8|. Likewise, it was pointed out by Barrow and 
Tipler, [5] that there are similar anthropic restrictions on the magnitude of any 
cosmological constant, A. If it is too large in magnitude it will either precipitate 
premature collapse back to high density (if A < 0) or prevent the gravitational 
condensation of any stars and galaxies (if A > 0). Thus, we can provide good 
anthropic reasons why we can expect to live in an old universe that is neither 
too far from flatness nor dominated by a much stronger cosmological constant 
than observed (|A| < 10|Aobs|)- Our results for varying a suggest that there 
might be significant anthropic constraints if A or the spatial curvature is too 
small to prevent domination before atomic structures become impossible. 

3 Observations in Space and in the Lab 

Studies of relativistic fine structure in the absorption lines of dust clouds around 
quasars by Webb et al., [21 [3], have led to widespread theoretical interest in 
the question of whether the fine structure constant has varied in time. The 
quasar data analysed in refs. [21 13] consists of three separate samples of Keck- 
Hires observations which combine to give a data set of 128 objects at redshifts 
0.5 < z < 3. The many-multiplet technique finds that their absorption spectra 
are consistent with a shift in the value of the fine structure constant between 
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these redshifts and the present of Aa/a = [a{z) — a]/a — —0.57 ± 0.10 x 
10~^, where a = a{0) is the present value of the fine structure constant [2 [3]. 
Extensive analysis has yet to find a selection effect that can explain the sense 
and magnitude of the relativistic line-shifts underpinning these deductions. A 
smaller study of 23 VLT-UVES absorption systems between 0.4 < z < 2.3 
by Chand et al. [M] initially found Aa/a = -0.6 ± 0.6 x 10"^ by using an 
approximate version of the full MM technique. However a recent reanalysis of 
the same data by Murphy et al. using the full unbiased MM method increased 
the uncertainties and suggested the revised figure of Aa/a = — 0.64±0.36 x 10~^ 
for the same data [3]. 

Any variation of a today can also be constrained by direct laboratory searches. 
These are performed by comparing clocks based on different atomic frequency 
standards over a period of months or years. Until very recently, the most strin- 
gent constraints on the temporal variation in a arose by combining measure- 
ments of the frequencies of Sr [35j, Hg-|- [36], Yb-I- [37], and H [38] relative to 
Caesium: a/a = (— 3.3±3.0) x 10~^^yr~^. Cingoz et al. also recently reported 
a less stringent limit of a/a = — (2.7± 2.6) x 10~^^yr~^ [3^; however, if the 
systematics can be fully understood, an ultimate sensitivity of lO^^^yr"^ is 
possible with their method [40j . If a linear variation in a is assumed then the 
Murphy et. al. quasar measurements equate to a/ a = (6.4 ± 1.4) x 10"^® yr^^. 
If the variation is due to a light scalar field described by a theory like that of 
BSBM [T3], then the rate of change in the constants is exponentially damped 
during the recent dark-energy-dominated era of accelerated expansion, and one 
typically predicts a/a = 1.1 ± 0.3 x 10^^^ yr^^ from the Murphy et al data, 
which is not ruled out by the atomic-clock constraints mentioned above. For 
comparison, the Oklo natural reactor constraints, which reflect the need for the 
Sm^^^ + n ^ Sm^^" + 7 neutron capture resonance at 97.3 meV to have been 
present 1.8 — 2 Gyr (z — 0.15) ago, as first pointed out by Shlyakhter [24], are 
currently Aa/a = (-0.8 ± 1.0) x 10"* or (8.8 ± 0.7) x 10"^ (because of the 
double-valued character of the neutron capture cross-section with reactor tem- 
perature) and Aa/a > 4.5 x 10~^ (617) when the non-thermal neutron spectrum 
is taken into account. However, there remain significant environmental uncer- 
tainties regarding the reactor's early history and the deductions of bounds on 
constants. The quoted Oklo constraints on a apply only when all other con- 
stants are held to be fixed. If the quark masses to vary relative to the QCD 
scale, the ability of Oklo to constrain variations in a is greatly reduced ^l]- 

Recently, Rosenband et al. |42j measured the ratio of aluminium and mer- 
cury single- ion optical clock frequencies, /Ai+//Hg+, repeated over a period 
of about a year. From these measurements, the linear rate of change in this 
ratio was found to be (—5.3 ± 7.9) x 10~'^^yr~^. These measurements pro- 
vides the strongest limit yet on any temporal drift in the value of a: a/a = 
(— 1.6± 2.3) X 10^^^ yr^^. This limit is strong enough to strongly rule out the- 
oretical explanations of the change in a reported by Webb et al. [S] [3] in terms 
of the slow variation of an effectively massless scalar field, even allowing for the 
damping by cosmological acceleration, unless there is a significant new physi- 
cal effect that slows the locally observed effects of changing a on cosmological 
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scales (for a detailed analysis of global-local coupling of variations in constants 
seeRefs. [4^). 

It has been noted that if the 'constants' such as a or ^ can vary, then 
in addition to a slow temporal drift one would also expect to see an annual 
modulation in their values. In many theories, the Sun perturbs the values of the 
constants by a factor roughly proportional to the Sun's Newtonian gravitational 
potential (the contribution from the Earth's gravitational potential is about 
14 times smaller than that of the Sun's at the Earth's surface). Hence the 
'constants' depend on the distance from the Sun. Since the Earth's orbit around 
the Sun has a small ellipticity, the distance, r, between the Earth and Sun 
fluctuates annually, reaching a maximum at aphelion around the beginning of 
July and a minimum at perihelion in early January. It was shown in Ref. [44j 
that in many varying constant models, the values of the constants measured 
here on Earth, would oscillate in a similar seasonal manner. Moreover, in many 
cases, this seasonal fluctuation is predicted to dominate over any linear temporal 
drift [g. 

Specifically, let us suppose that the Sun creates a distance-dependent per- 
turbation to the measured value of a coupling constant, C, of amplitude 6hiC = 
C{r). If this coupling constant is measured on the surface of another body (e.g. 
the Earth) which orbits the first body along an elliptical path with semi-major 
axis a, period Tp, and eccentricity e ^ 1, then to leading order in e, the annual 
fluctuation in C, (JCannuai will be given by 

= -cccos j +0(e ), (10) 

where cc = eaC'{a), C'{a) = dC{r) / dr\r=a and t = nTp, for any integer n, 
corresponds to the moment of closest approach (perihelion). In the case of the 
Earth moving around the Sun, over a period of 6 months from perihelion to aphe- 
lion one would therefore measure a change in the constant C equal to 2cc ■ Using 
5 ln(/Ai+//Hg-i-) = (3.19-)-0.008)fo/a, [i^ . a maximum likelihood fit to the data 
gives Ca = ea6a'{a) = (-0.89 ±0.84) x lO'^^, where a = 149, 597, 887.5 km is 
the semi-major axis of the Earth's orbit, and Sa{r) is the perturbation in a due 
to the Sun's gravitational field. Assume that over solar system scales, the val- 
ues of the scalar fields on which values of the 'constants' depend, vary with the 
local gravitational potential. Hence, we have Sa{r)/a = kaAUQ{r), where ka is 
a theory-dependent multiplier, AUq is the change in the gravitational potential 
of the Sun: [/©(r) = -GMq/t, and so eaAC/^(a) = eGMQ/a = 1.65 x 10^1°. 
Hence, we find: 

fc^ = (-5.4±5.1) X 10^^ (11) 

The frequency shifts measured by Rosenband et al. [42] were not sensitive 
to changes in the electron-proton mass ratio: fi = me/wp. Measurements of 
optical transition frequencies relative to Cs, Refs. [351 [SZl [35] j are sensitive 
to both fj, and a. H-maser atomic clocks [45] are also sensitive to variations 
in the light quark to proton mass ratio: q — mq/nip. We can use all these 
observations if we define two more gravitational coupling multipliers, fc^ and fcq. 
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by 5pL/ii = kf,AUQ, and Sq/q = fcqAJJ©. Refs. [351 [Ml SS] give fc„ + 0.36fc^ = 
(-2.1 ± 3.2) X 10^6^ ko, + O.Uk^ = (3.5 ± 6.0) x 10"^ and k^ + O.Uk^ = 
(1 ± 17) X 10~^ respectively. We also performed a bootstrap seasonal fluctuation 
fit (with 10^ resamplings) to the Yb^ frequency measurements of Peik et al. 
[57] giving ka + 0.51fc^ = (7.1 ± 3.4) x 10^*^ [?]. Combining these bounds with 
Eq. (HI]) gives /c^ = (3.9 ± 3.1) x 10-^ fcq = (0.1 ± 1.4) x 10"^. Recently, 
Blatt et al. [35] combined data from measurements of H-maser [45] and optical 
atomic clocks [35] [33], to bound the multipliers, ka, k^ and kq, finding ka = 
(2.5 ± 3.1) X 10~^.The constraint on ka derived in this paper from the data 
of Rosenband et al. [42] therefore represents an improvement by about two 
orders of magnitude over the previous best bound. This improved bound on 
ka combined with data found by Peik et al. [37] has also produced an order of 
magnitude improvement in the determination of A;^ and a slight improvement 
in the constraint on kq. 

Seasonal fluctuations are predicted by a varying constant theory because 
the scalar field which drives the variation in the constant couples to normal 
matter. The presence of the Sun therefore induces gradients in scalar fields 
and associated varying 'constants', and it is essentially these gradients that are 
detectable as seasonal variables. As mentioned earlier, gradients in a scalar 
field which couples to normal matter result in new or 'fifth' forces with pseudo- 
gravitational effects. In the case of varying a and ^ theories, these forces are 
almost always composition dependent, which would violate the universality of 
free-fall and hence the weak equivalence principle (WEP). The magnitude of 
any composition-dependent fifth force toward the Sun is currently constrained 
to be no stronger than 10~^^ — 10""'^^ times than the gravitational force [46]. In 
the context of a given theory the constraints from WEP tests indirectly bound 
ka- Indeed, they often provide the tightest constraints on ka [?. IT7lH5] . 

A recent thorough analysis of the WEP violation constraints on ka [47] |49] 
found ka = (0.3± 1.7) x 10^^, with a similar constraint on kq. It must be noted, 
however, that this result is still subject to theoretical uncertainty, especially 
regarding the dependence of nuclear properties on quark masses. For instance, 
it was also noted in Ref. [17] that if certain (fairly reasonable) assumptions 
about nuclear structure are dropped, the la error bars on ka increase by about 
an order of magnitude to: ±1.4 x 10~*. Despite these uncertainties, for many 
theories of varying a, WEP violation constraints from laboratory experiments 
or lunar laser ranging [50] still provide the strongest, albeit indirect, bound on 

ka- 

4 Conclusions 

We have described some of the history of the study of varying constants in 
physics. This area of research has been reinvigorated by a significant body of 
observational data, drawn from quasar absorption spectra, which is consistent 
with a change in the value of the fine structure constant, a, over of a few parts 
in a million over 10 billion years. So far, these data have neither been reliably 
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confirmed nor contradicted by other observational studies and this confrontation 
is keenly awaited. We described how a simple self-consistent theory of varying 
a developed by Sandvik, Barrow and Magueijo can be constructed and the clear 
pattern of variation that it predicts in the universe: no variation of a during 
the radiation era and a logarithmic time increase during the cold-dark-matter 
era, followed by a resumption of no time variation in a after the universe begins 
to accelerate during the dark-energy era. A fit of these simple models to the 
observational data fixes the one free parameter defining the theory and pre- 
dicts a violation of the weak equivalence principle at a 10~^^ level that is easily 
detectable from space. We also described exciting new developments in the lab- 
oratory search for varying a. These experiments are for the first time achieving 
the sensitivities of the indirect astronomical bounds and in the next few years 
we may be able to draw some strong conchisions about varying constants from 
the confluence of laboratory experiments and large new astronomical data sets. 

Acknov^rledgements I am very grateful to Mariusz Dabrowski for his help 
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